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Large Chern number phases in Haldane model become possible if there is a multiplication of 
Dirac points in the underlying graphene model. This is realized by considering long distance hop- 
ping integrals. Through variation of these integrals, it is possible to arrive at super-merging band 
touchings, that up to N7 graphene are unique in parameter space. They result from the synchro- 
nized motion of all supplementary Dirac points into the regular ±K points of graphene. The energy 
dispersion power law is usually larger then the topological charge associated to them. Moreover, 
adding distant-neighbor hoppings in the Haldane mass allows to sweep large number Chern phases 
in the topological insulator. 



I. INTRODUCTION 

The Haldane modeP^ is the first topological insulator 
that presents quantum Hall effect without an external 
magnetic field. It is a two-band system with bands that 
has a nontrivial topology. The bands are characterized by 
Chern numbers that are proportional to the conductance 
carried by edge states. This model was a playground for 
ideas that eventually led to the pre diction and discovery 
of the Z2 topological insulators .^^El 

Here we revisit the Haldane model and show in prac- 
tice how the addition of hopping integrals between dis- 
tant neighbors can lead to a multiplication of topological 
phases with a large Chern number. Recent work suggests 
a way to produce a flat bands with arbitrary Chern num- 
bers in multi-layer systems.'^^ In contrast, we constrain 
ourselves to the two-bands Haldane system and we do not 
seek flatness of bands. Admittedly this is not a very phys- 
ical way to increase the topological index characterizing 
a band, because the contribution from distant neighbors 
are small. The conceptual advantage is that one can fully 
describe the phase diagram of such systems and analyti- 
cally predict its topological transitions. The system can 
be understood from a decomposition of the model in an 
underlying gapless graphene model and a Haldane mass. 
The variation of the topological index is related at once 
to the multiplication of nodes in the energy dispersion 
for the gapless system and to the rapid oscillations in the 
Haldane mass term. In general, if the two-band under- 
lying gapless system admits 2n Dirac points, then the 
Chern number can vary from —n to n'.^ 

In Sec. [H] we present the theoretic tool to compute 
analytically the Chern number in a two-band system. 
That allows to immediately discriminate the topological 
phases. In Sec.|IV[ we treat the underlying graphene with 
long-distance hopping integrals. Up to N7 (next x 6 near- 
est neighbor) graphene, we investigate the multiplication 
of Dirac points through the addition of long range hop- 
ping. We also show that there are unique super-merging 
points where the all Dirac points merge to ±K points in 
the graphene Brillouin zone (BZ). Satellite Dirac points 
can be found by perturbing around these special band 



touchings. The topological charge associated to them 
can be immediately established from a sum over the ad- 



ditional Dirac points. In Sec. IV we consider the effects 
of gapping the Dirac points with a Haldane mass term. 
The phase diagram for the modified Haldane model are 
shown. 



II. CHERN NUMBER IN TWO-BANDS 
MODELS 



Generic two bands Hamiltonian on a two dimensional 
Bravais lattice reads 

^=^1/'^ E /^«Mk)clkC/3k (1) 



a, (1=1,2 



with c|j]^ the creation operator of Bloch state with wave 
vector k and where a constitutes a pseudospin index re- 
sulting from either two sublattices or two orbitals per 
unit cell. The elements ha/iOi) form a 2 x 2 Hermitian 
matrix h(k.) that can be written 



h(k) 



^/l^(k)cr^. 



(2) 



with ctq the identity matrix and cri^2,3 the Pauli matrices. 
/ip=o,3(k) comes from intra- sublattice contributions haa 
and h^^i 2(k) from inter- sublattice contributions h^p- 
The real valued functions /i/j(k) can be further split into 
even and odd components /i^(k) — h'^{'k) + /i^(k), where 



Kik) = Ki-k) and K{k) 



-h°{~k). For time- 



reversal symmetric, spinless Hamiltonians, ^;i=o,i,3(k) 
are purely even and /i2(k) purely odd. 

The spectral decomposition of matrix h{k) reads 



(3) 



with band energies e±(k) = ho{k) ± |h(k)| and eigen- 
band projector P±(k) = ^{ao±h- <T/|h|) where h(fc) — 
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(/ii, /i2, ^s)- Component /io(k) breaks particle- hole sym- 
metry by shifting the energy bands, but does not inter- 
vene in the gap |h| or in the projectors P± which deter- 
mine the topological properties of the Hamiltonian. 

In an insulating phase the three components of h(k) 
never vanish simultaneously and |h| remains finite for 
any k. Such insulating phase can be characterized by a 
topological index, the first Chern number C, that shows 
how many times the Brillouin zone wraps around the unit 
sphere traced by h = h/|h|. One can choose to index the 
system with the Chern number associated to the lowest 
band e-(k) 

C=-^ / d^kh-{dkhxdkh), (4) 

where the integral is over the Brillouin zone. The Chern 
number is zero unless one of the component of h(k) 
breaks time-reversal symmetry. Furthermore, a non-zero 
value of C requires that any submodel obtained by con- 
sidering only two components of h must exhibits band 
touchings at some finite set of points in the BZ.^^In fact, 
an explicit calculation of C is made easy by considering 
such a gapless system containing only two components 
of h. When the band touchings of the gapless submodel 
have linear dispersion (i.e. they are Dirac points), C can 
be calculated by treating separately the chirality of the 
Dirac points and the sign the third component of h (the 
mass term that gaps the system) at these Dirac points. 
The Chern number then read^ 

c = ^ E (5) 

where Di are the set of Dirac points for a simplified two 
components model without hi term, and 

Vk e A, X*(k) = sgn{dkh X dk h)^ (6) 

are their corresponding chirality. Such formula permits 
an analysis of the Haldane model by separately study- 
ing the underlying gapless graphene model and the sign 
of the mass term. The caveat of the above formula is 
that it works only for point band touchings with linear 
energy dispersion (Dirac points). However we shall see 
that there can be as well Fermi lines or point band touch- 
ings with higher power dispersion. The latter ones can be 
understood as merging of many Dirac points. Then the 
topological charge of the merging points is just the sum 
of the chiralities for the Dirac points that are converging 
to it. This fast calculation of charge associated to a band 
touching will be referred to as the sum rule.^^ 

III. DISTANT-NEIGHBOR HOPPING IN 
GRAPHENE 

As seen in previous section, the possibility of a two- 
band insulator with a Chern number C = n requires to 




FIG. 1. The possible hoppings in graphene N9 model. From 
a central B atom, the neighbors are arranged in concentric 
circles. The hopping integrals from the central B atom to a 
site placed on a circle is denoted by n growing with the 
distance between sites. 

build a reduced two-band gapless model with at least 2n 
Dirac points. Let us consider the extended tight-binding 
graphene model, including distant Nn (nextx(n — 1)- 
nearest-neighbor) hopping terms on the hexagonal lat- 
tice. The eventual Dirac point will eventually be gapped 
by a Haldane mass to yield a topological insulator with 
large Chern numbers. The hexagonal lattice is a bipar- 
tite lattice built out of two interpenetrating triangular 
Bravais sublattices A and B. Let us denote by t„ the 
(isotropic) Nn the intra- and inter-sublattices hopping 
terms (see Fig.[l]). The parameter ti corresponds to usual 
nearest-neighbor Nl graphene. In this section we consid- 
ered only inter-sublattice hoppings t„ in units of ti , such 
that there are n — 1 free parameters. The intra-sublattice 
hopping terms contribute to the identity Pauli matrix ctq 
and are neglected in the following. For real- valued the 
matrix /i(k) preserves time- reversal and inversion sym- 
metries (e.g. /i*(— k) = h(k) and o'i/i(— k)^! — /i(k)). 
Moreover, when considering only inter-sublattices hop- 
ping i„, there is a sublattice symmetry characterizing 
the system. The symmetry is represented the operator 
CT3 which anticommutes with /i(k).^ Explicitly, /i(k) reads 

^('^) = (r?k) ^0' 

with /(k) = /ii(k) — i/i2(k) or 

/(k) = E*»fn(k), (8) 

n 
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where the functions 5n(k) up to N9 are tabulated in Tab|l] 
in which ai = -\/3a(^, ^) and a2 = y/3a[ — 5, ^) de- 
note the primitive vectors of the triangular sublattices. 
The hexagonal lattice constant a is set to 1 from now 
on. The explicit form of the function ,9ri(k) corresponds 
to a Bloch basis such that (7n(k + G) = 5,1 (k) with G a 
reciprocal lattice vector. 



A. Nearest-Neighbor Nl graphene review 

Before studying Nn graphene, let us briefly review the 
usual properties of Nl graphene where /(k) = tigi(k). 
The two energy eigenvalues are given by e±(k) — ±|/(k)| 
and there is a gap separating the two bands. Band 
touchings occur at isolated positions ±K correspond- 
ing to zeroes of /(k). For Nl graphene the zeroes cor- 
respond to the two nonequivalent Brillouin zone cor- 
ners ±K = ±(aj - a*)/3 where a* = |^(^,|) and 

a2 = |^(— ^,|). At each of these points, there are 
two degenerate zero-energy eigenstates. As illustrated 
in Fig. [2| the bipartite property allows to project one 
state entirely on the A sublattice and the other on the B 
sublattice.'^ Altogether there are four zero-energy states, 
each labeled by two indices; a valley index correspond- 
ing to ±K and sublattice index A or B equivalent to the 
eigenvalues ± associated to sublattice symmetry opera- 
tor CT3. The time-reversal transformation exchanges val- 
ley index without changing sublattice index. Inversion 
(represented by cti) exchanges both valley and sublat- 
tice indices. Hence the product of the two operations 
exchanges sublattice index only. 

In the neighborhood of the band-touchings ±K, one 
can expand the function /(k) in small momenta q = 
q'(cos6i,sin0). It follows that /(±K -h q) ~ ge'F'^ and 
the linearity in q identifies the band-touchings as mass- 
less Dirac fermions. More generally, if the band touch- 
ing has / oc ((/e^*^)" then its respective chirality is n. 
This translates the fact that the two dimensional vector 
(/ii,/i2) oc {cos{n6) , sm{n9)) rotates counter-clockwise 
by 27rn for 6 sweeping once the interval [0, 27r). Here, 
for n = 1, it follows immediately that Eq. ^ and the 
low-energy expansion, both predict x(±K) = ±1. 



B. Band touchings in Nn graphene 

Let us study how are the properties of the zero-energy 
states modified when distant neighbor hoppings is consid- 
ered. For Nn graphene, band touchings occur at positions 
±k corresponding to zeroes of /(k) = J^n^ngnO^)- Pre- 
vious solutions, ±K, obey (7„(±K) = and thus remain 
zeroes of /(k) regardless of the new hopping integrals 
tn>i- To find the positions of other zeroes, one can keep 
k on the three high-symmetry lines T joining F, ±K, and 
M. These lines are globally invariant under time reversal, 
C3, C2 and inversion with respect to F point. Without 




FIG. 2. Real space representation of the four zero-energy 
eigenstates of Nl graphene.® Filled (empty) bullets repre- 
sent A (B) sublattice sites. Wave functions components 
on different lattice positions are related by Bloch theorem 
?/'k(r + R) = e''''^'(/)k(r) with R any Bravais lattice vector. 
Let us denote e = e'^^^'^i = e^'^'/^ then e^'^"'^ = t with 
1 -f e + e* =0 and — — \ Left figure corresponds to 
valley K and A sublattice and right figure to valley K and 
B sublattice. The wave functions in — K valley are obtained 
by complex conjugating the amplitudes at K. Wave function 
amplitudes are invariant under C3 rotation around a lattice 
site and under translations (Rx = m(ai + a2)) perpendicu- 
lar to K and exhibit periodicity under translations parallel to 
K with a period Ry = 3(ai — a2). Time reversal exchanges 
valley index without changing sublattice index. Inversion ex- 
changes both valley and sublattice indices. The product of 
the two operations thus exchange sublattice index only. 

loss of generality, let us analyze the T line, k = fc(l,0) 
(see. Fig. [4|. Along it 5n(k) is a real function and the 
condition /(k) ~ translates into a polynomial equation 

hi{x) — for the variable x — cos(^fc). Up to N8 this 
polynomial reads as 

hi{x) = A(yX + -^p{x), with 

p{x) = 4(t7 -I- tf,)x^ + 2{U - h)x^ + {h - 4^7 - ?,ts)x 
+ 2-f-i4 + t7 + f (9) 

The maximum number of solutions is given by the de- 
gree of the polynomial and it clearly increases with the 
range of hopping, but not systematically since N7 and 
N8 correspond to a polynomial with the same degree. 
When all the solutions are distinct, they correspond to 
band touchings with linear dispersion in kx direction. As 
anticipated x = — ^ (at ±K) is a solution regardless the 
value of tn . Other physically meaningful solutions must 
verify |a;| < 1. For each such solution upon applying a C3 
rotation, one can associated three band touching points 

at ki = fc(l,0), ka = k{-\, ^) and kg = k{-\,-^). 
From time-reversal symmetry, it follows that there are 
three additional touching points at — ki_2.3- Hence a so- 
lution \x\ < 1 with X 7^ {±1,-1/2} implies at least six 
nonequivalent band touching points at ±ki.2,3. In con- 
trast, a solution x = — 1 is associated to a single touching 
point at F, a solution a: = — ^ to the two nonequiva- 
lent BZ corner ±K and a solution a; = 1 to the three 
nonequivalent M points. In a Nrt graphene model with a 
polynomial hi (x) of degree m < n the maximum number 
of nonequivalent band touchings points per valley on a T 



4 



, . physical chemical . , , ^ 

Nn hoppmg ^.^^^^^^ ^.^^^^^^ neighbors s„(k) 

In Tl i i 3 l + g-ik.ai ^g-ik-aa 

N3 t3 2 3 3 -"^2) +e*-{a2-ai) _,_g-ik.(ai+a2) 

N4 ^ 3 g e*''''*! + e*''''^2 _(_ e~2»k-ai _|_ g-2ik-a2 _|_ gik'(ai-2a2) ^ gik-(a2-2ai) 

N7 y/l3 5 g gik-(2ai-a2) _j_ gik-(2a2-ai) ^ g2ik-{ai-a2) _|_ g2ik'(a2-ai) ^ g-ik-(2ai+a2) _|_ g-!k-(2a2+ai) 

N8 tg 4 5 3 gik-(ai+a2) _|_ gik.(ai-3a2) ^ gik-{a2-3ai) 

N9 5 g g-3ik-ai _|_ g-3ik a2 ^ gik-(2ai-3a2) _|_ gik-(2a2-3ai) _|_ g2ik ai ^ g2ik-a2 



TABLE I. Properties of Nn AB inter-sublattices hopping terms. Physical distance is counted in units of lattice constant a. 
Chemical distance is the smallest number of bonds passed while hopping between two sites. In the "neighbors" column are 
counted the number of sites at a given physical distance from a chosen central site. In contrast, note that any site has 3n 
neighbors located at a chemical distance n from it. The primitive vectors of the triangular sublattice are ai = V3a{l, ^) and 
a2 = V3a(-l,^). 



lines is thus (1 + 3(m — 1)). For example, in N7 graphene 
the degree of the polynomial is m = 4, and therefore 
there are a maximum of 10 Dirac points per valley. As 
a final remark concerning the band touching points, we 
stress that we cannot exclude the possibility of having 
additional touching points outside the high-symmetry T 
lines. However, the only case encountered in the numer- 
ical simulations is that of Fermi lines (zero energy lines) 
which connect the zero energy solutions located on the 
T lines. These are particular solutions that can be ex- 
pected when a non-degenerate zero on the T line exhibits 
a vanishing chirality (see Sec. IIIF for an example). 



C. Zero-energy state wave functions in Nn 
graphene 

Similarly to Nl graphene, there are two degenerate 
zero-energy eigenstates that correspond to each band- 
touching point of Nrt graphene. The bipartite property 
is still valid for Nn graphene and it allows to project one 
zero-energy state on the A sublattice and the other on 
the B sublattice. The real space representation of these 
two energy states is illustrated in Fig. |3] for a generic 
k = fc(l,0) on a T line. The wave function exhibits 
translation invariancc in the direction perpendicular to 
k and it is multiplied by a phase — e^^^'^^ on both A 
and B sublattices when translated by one unit along k. 
Fig. [3] is especially useful as it allows to quickly construct 
the polynomials hi{x) at all orders. 




FIG. 3. Real space representation of a generic zero-energy 
eigenstates of Nn graphene projected on the sublattice A (in 
valley K). Wave functions components on different lattice po- 
sitions are related by Bloch theorem i/jk(r -I- R) = e^^'^ipk{r) 
with R any Bravais lattice vector. For k on high-symmetry 
lines we have z = e'"""! and z* = 6*""^ (for k / ±K, 
1 + z + z* 7^ and 2^ 7^ 1 7^ z*'^). The three additional states 
are obtained by performing a C3 rotation around the center 
of a hexagon. Solutions in opposite valley are obatained by 
complex conjugating the amplitudes. A similar picture can 
be drawn for states projected of the B sublattice. 



D. Velocities and chirality of Dirac points in Nn 
graphene 

When all the band touching points on a given T line 
are distinct, each one of them may correspond to a Dirac 
point k. The energy dispersion near the touching point 
is obtained by expanding to first order in small momenta 
q = g(cos 9, sin 6) 

f{k + q)^qidkMcose-idk^h2Sm9) + Oiq'), (10) 



where on the T line (k = (fc,0)), we used the property 
dk hi — dk^h2 = 0. Let us define the velocities Cx — 
dk^hi and Cy — dkyh2- A band touching point k is a 
Dirac point, if both velocities are non- vanishing at k, 
7^ 7^ Cy. More quantitatively c^, Cy read as: 

cx(x) = T2\/3Vl-a;2[p(a;) + {x + \)p' {x)], (11) 
cy{x) = -i[p{x) - [x +\)r{x% (12) 
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where the sign =F of Cx corresponds to band touchings 
associated to ±K valley, and the polynomial r{x) is given 
by 



r{x) = 16^82;^ + 4(t4 - ts)x'^ + 4(i4 - 4ts)x 
+l + t3-6t4 + 2tj + 5ts. 



(13) 



To simplify the above equations, it is opportune to study 
separately the velocities for band touchings at the ±K 
points (x = —1/2), and eventual solutions away from the 
±K points on the T line for x ^ ~l/2 and p(x) = 0. 
Let us take the band touchings only at the K valley, 
knowing that the changes sign at the opposite valley. 
The corresponding velocities read 



--(l-2t3-t4 



5*7 



4*8), 



(14) 



-2Vl5VT 
3(a: 



.x2(.T+ i)p'(a;). 



)r{x). 



The above equations indicate that the two velocities are 
equal in magnitude and eventual Dirac points will have 
isotropic cones at ±K. Also note that a merging of Dirac 
points in K valley creates a energy dispersion of higher 
order in q and these is equivalent to vanishing of the 
velocities to Cx^y\K_ = 0. At band touchings different 
from ±K one can use the condition p(x) — to simplify 
the expression of the r(x) polynomial: 

r{x) = 2(l-x)[4(*7-i8)a;^+4(i7-*8)2: + *3-2t4+i8]- 

(15) 

At time-reversal points, Y[x = 1) and M(x = —1), the 
velocity Cx is always zero. F point (center of the BZ) is 
the band bottom and presents an isotropic energy disper- 
sion; therefore Cy vanishes together with Cx (as seen from 
Eq. p^). In contrast, at the M point, Cy is not nec- 
essarily zero. For example, this allows in N3 graphene 
for M band touchings with linear dispersion in fc^ and 
quadratic in ky. These semi-Dirac points correspond to 
a merging of two Dirac points with opposite chirality. 

If all the band touchings are Dirac points, then their 
chirality Q follows from Eqs. 14 



x + 



-2 :x(±K) 



±1, 



Tsgn[p'(a;)r(a::)], (16) 



where ztk^ denote the position of the additional Dirac 
points associated to the ±K valley. The next sections 
exemplifies the above theory to the concrete cases of N3 
and N4 graphene. 



satellite band touching points orbiting around each reg- 
ular Dirac point (±K). Indeed solving Eq. (|9| (with 
^4 = ^7 = ^8 = 0), it follows that in addition to solution 
X — (at ±K), there is a solution x — which may 
give rise to up to six touching points at 



±ki = ±fc(l,0), 

1 



±k; 



2' 2 



±k2 = ±k 



1 ^/i 
2' ~Y 

2fa 



(17) 



where ik^ points are associated to the ±K valley. Phys- 
ically meaningful solution corresponds to < 1 and 
hence have an existence domain given by 



is e (-00, -1) U (1/3, CX)). 



(18) 



For t^ e (— 1,— oo) satellite touching points appear at 
F (t3 ~ —1) and move along T line and reach the E 
point (x = 1/2, <3 = — cxd), midway between K and F. 
For ^3 e (1/3, oo) satellite touching points appear at M 
(t3 — 1/3) and move along the T line and reach again 
S (x = 1/2, = oo) (see Fig. |4|. For ^3 ^ 1/2, the 
satellites are Dirac points away from the regular Dirac 
points ±K, X — 7^ ^1/2. The chirality associated 
to the three satellite Dirac point ki^2.3 in valley K reads 



1 



2t. 



-Sgn(i3(l + t3))- 



(19) 



The chirality x is always opposite to points associated to 
the — K valley. 

As already emphasized,ESl there is a particular value, 
^3 = 1/2, that corresponds to a merging of three satel- 
lite Dirac points with a central regular Dirac point. This 
case is realized when x = —1/2 is a double root of the 
polynomial hi{x). Here p{x) — and therefore the ve- 
locities Cx^y vanish simultaneously, indicating the for- 
mation of a band touching with a higher than linear 
dispersion. Note however that at the merging point, 
x{x — —1/2) = sgn{cxCy) is not well-defined. Never- 
theless, from Eq. ( 19 ) it is apparent that the satellite 



points close to merging at ±K have an opposite chirality 
from the central Dirac point x(iK) = ±1. Then the sum 
rule dictates that the chirality at the merging point is the 
sum of chiralities over the colliding Dirac points. At±K- 
merging this yields: xih = 1/2) = ±(1 - 3) = =f2.I^ 

The chirality of the merging point can be equally de- 
termined by expanding the energy dispersion at ±K. It 
suffices to find it at K, knowing that time-reversal sym- 
metry demands opposite chirality at — K. Expanding at 
^3 = 1/2 in small momenta q — q{cos6,sm6) it follows 
that 



E. Dirac points and merging for N3 graphene 

The isotropic N3 graphene was already investigated 
in Ref. [ini Here the presence of a sufficiently strong 
^3 hopping integral was shown to produce three more 



/(±K + q) 
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(20) 



This indicates that the band touching at the merging of 
all the Dirac points in a valley has a quadratic dispersion 
and a topological charge of ^2 in valley ±K. 
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(b) 
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FIG. 4. (Color online.) Evolution in BZ of a satellite Dirac 
point in N3 graphene on the high-symmetry T-line: T—K.—M. 
The evolution of the satellite point is represented in blue when 
t3 varies from — oo to —1 and in red, when ts varies from 1/3 
to 00. 



F. Dirac points for N4 graphene 

For N4 graphene, solving Eq. Q yields, beside the 
solution X — —-^ (at ±K), two additional solutions x± = 
T^[t3 ± {tl + 8tl + AtiU - 4t4)^/^l (x^ < x^) such 



^ yo^ -r 01.4 -r 4i4t3 - Ati)'-''^] (x^ 

that there are up to 7 band-touching points per valley. 
More quantitatively, for < t^,t4 < 1, one obtains the 
existence domains for additional solutions when |a;± | < 1. 
Explicitly, |a:+| < 1 for 



4 > — and 2( 
^ - 10 



(-4 



t4) < <3 < 



1 -I- 2tA 



(21) 



Similarly, < 1 holds for 



1 , l + 2U\ 

Ia < — and ti, > or 

10 3 J 



(22) 



The existence domains are represented graphically in 
Fig. [5] Note that the two solutions coexist when ^3 < 
i^i^P*. In the coexistence region one can expect gener- 
ally to have 7 Dirac points per valley. 

In their existence domain, Eq. ( 16 1 determines the chi- 
rality in the K- valley 



X{x±) = ±sgn(t3 - 2ti). 



(23) 



However the chirality information is exact when the so- 
lutions x± stand for Dirac points. The model presents 
a rich phenomenology and the investigation of the so- 
lutions indicates that for particular parameters there are 
also band touchings different from the simple Dirac point 
case. 

Remember that each solution x± stands for a triplet 
of solutions at each valley. Then there are different sce- 
narios for the behavior of Dirac points. There are cases 
similar to the N3 graphene where there is a single triplet 
of solutions merging to the central Dirac point to yield a 
point with high-energy dispersion. There are cases where 



FIG. 5. (Color onhne.) (a) < 1. (b) |a::+| < 1. Existence 
domains and corresponding chirality of solutions |a;±| < 1 in 
{ti, ta) parameter space and in the K valley. The region with 
positive (negative) chirality is represented in red (blue). The 
green line ts = 2t4 where chirality changes is associated to the 
existence of Fermi lines instead of Dirac points. The super- 
merging point ts — and t4 = | at the intersection of ts = 2t4 
line with the domain border curve t-j = 2{^t4 — t\ — 14) is 
indicated in yellow. 



the two triplets merge with one another to yield a new 
triplet of band touchings with quadratic dispersion in one 
direction and linear in the other. There is also a unique 
super-merging point where all Dirac points in a valley 
merge. A completely new feature to the phenomenology 
of band touchings in N4 graphene is the formation of 
Fermi lines for specific values of parameters. 

The first case is that of a line in parameter space where 
only one triplet given by the |a;± | -solutions merges with 
the central Dirac point. These are obtained under the 
condition that either x+ = ^1/2 or x^ = —1/2, 



x± 



1 - ti 



(24) 



A different case is that of the triplet satellite Dirac 
points merging two by two to form semi-Dirac points, 
i.e. band touchings with quadratic dispersion in the di- 
rection of merging and linear in the direction perpendic- 
ular to it .133 They correspond to a scenario where two 
Dirac points with opposite chirality collide. From the 
condition a;+ ~ x- , they are determined on the line 
is = 2{^Jt4 — t\ — ti). This case is represented in Fig. [g] 

At the intersection of line ^3 — 21^ with the domain 
border curve ^3 = 2{yjt4{\ — ^4) — ti) {t^ — | and 
^4 — |), there is a super- merging point where there is 
a unique band touching per valley that can be under- 
stood as a collision of all additional Dirac points into 
the central (±K) one. This point in parameter space 
has a topological charge given by the sum of all Dirac 
point chiralities. Because the cancellation of the triplet 
charges the, final point will have a charge ±1 in the valley 
±K. Expanding in small momentum q — (/(cos 0, sin 0) 
around the super-merging point at ±K yields an effective 
/ function in K valley 



/(±K + g) 



27 

'40 



(25) 
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FIG. 6. Cross section slice through the energy dispersion of 
the conduction band near zero energy. A triplet of semi-Dirac 
points is formed around the cen tral Dirac point at K. On the 
parameter line = 2{^t4 — tl — ti) (here with — 0.25) the 
triplets of satellite points merge on the high-symmetry lines 
to form the semi-Dirac points. 



This result reinforces the sum rule calculation by showing 
a band-touching with cubic dispersion, but with a low 
topological charge ±1 at ±K. 

Finally, there is a phenomenologically new situation 
that is absent in the previously studied N3 graphene. 
Note that on the line — 2t4 the chirality (23) is zero 
even though there are nondegenerate solutions, x- =/= x^, 
away from the super- merging. This case corresponds to 
the existence of closed Fermi lines in the Brillouin zone 
that link the two solutions. Hence, in contrast with the 
cases studied until now, here the energy dispersion ex- 
hibits a line of zeroes outside the T line. One of the 
cases is represented in Fig. [7j where the a;± -solutions are 
connected by a Fermi line. However, even for a vanishing 
x+-solution, Fermi lines subsist and link band touchings 
associated only to a;_-solution. Aside from these numer- 
ical observations of the Fermi line at ^3 = 2<4, it remains 
a daunting task to analytically solve for general solutions 
away from the high-symmetry lines. However one can in- 
vestigate analytically the peculiarity of this case by con- 
sidering the behavior of the energy dispersion near the T 
lines. The absolute value of the energy for — 2ti is 



E = |4t4X^ + 4^42;?; - 4*4 + 1| V'4a;2 + 4xy + 1, (26) 

where x — cos{\/ikx/2) and y — cos{3ky/2). It is im- 
mediate to verify that the derivatives in the fcy-direction 
for a zero-energy solution x± on the T line /c(l, 0) vanish 
at all orders. This verifies that the solution x± are not 
longer point-like band touchings, but extend as Fermi 
lines in /c^-direction. 



-o-i^o.io 
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FIG. 7. An expansion in small momentum q around the K 
point of N4 graphene illustrates the formation of Fermi lines 
(lines of zeroes for the energy dispersion) around the regular 
Dirac point K in graphene, on the parameter line ts = 2*4. 
The hopping parameters are chosen near the super-merging 
at ts = 2/5 -I- 2(5 and t4 ^1/5 + 5 with 5 = 0.001. 



G. Super- merging at ±K in Nn graphene 

In the two preceding sections it was shown that for N3 
and N4 graphene it is possible to adjust the parameters 
^3,^4 so that for each valley all the additional touching 
points merge with the usual Dirac points at ±K (a super- 
merging point). This means that x = —1/2 is a double 
(respectively triple) zero of hi for N3 graphene, ^3 = 1/2 
(N4 graphene, ^3 = 2/5, *4 = 1/5). The possibility of 
finding a set of parameters *„ for which all the additional 
touching points merge with the usual Dirac points at K 
appears to be valid for all Nn graphene and relies essen- 
tially on the fact that the polynomial hi{x) is of a degree 
equal or less that the number n — 1 of free parameters i„. 
(More precisely, it can be proven that a model with hop- 
ping terms at a chemical distance m will result in poly- 
nomial hi{x) of maximum order to.) Note that because 
the number of free parameters grows faster than the de- 
gree of the polynomial there are no longer unique super- 
merging points for graphene Nm, with m > 7. At this 
super-merging the components Cx^y vanish and therefore 
one needs to go beyond a linear expansion to characterize 
the neighborhood of K. As an example, it was shown in 
Ref. JO, that for N3 graphene at the super-merging one 
obtains /(±K-|-q) ~ q^e^"^^^ which now identifies a gap- 
less quadratic dispersion, with a phase that is understood 
as resulting from the sum of the respective chirality of 
all the merging Dirac points. Similarly, for N4 graphene. 



it followed that /(±K -I- q) 



q^e 



=F*^ The location of 



the unique super-merging band-touching and their asso- 
ciated topological charge are given in Tab. [H] Note that 
the energy dispersion of super-merging band-touchings 
have a higher than linear dispersion. However the topo- 
logical charge of the converging triplets of satellite points 
is alternating and hence the resulting topological charge 
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FIG. 8. (Color online). The zero lines of /ii(k) = (in green) 
and ft2(k) = (in red) for N4 graphene. A small perturbation 
(+0.001) of t4 at the merging point = 2/5, t4 = 1/5 creates 
6 Dirac points around the stable Dirac point K. In the inset 
there is a zoom around K. The Dirac points are represented 
by full circles, •; there is a central K Dirac point in black, 
and 2 sets of satellite Dirac points, in blue and red. 

remains low. 

Finally, note that the above scenario of a unique super- 
merging together with alternating ±1 and +2 topologi- 
cal charge at ±K (see Tab. is not generally valid in 
Nn graphene, and in fact it already breaks down in the 
N8 model. For N8 graphene, the super-merging is no 
longer unique, but becomes a line in (is, ^4, ty, tg) pa- 
rameter space. Nevertheless, Eq. |9] implies that in N7 
and N8 graphene there is the same number of satellite 
Dirac points per valley, because p{x) remains a third or- 
der polynomial. An expansion near the super-merging 
point for N7 graphene (see Tab. |ll]) yields in K- valley 

97 

/(K + q) = -n* [n' - Utsin' - n*'))] , (27) 

with t: = Qx + iqy For vanishing t^, one recovers a topo- 
logical charge —2 for the band touching at K, at the 
super-merging in N7 graphene. However when reaches 
the critical value 1/12 the band touching clearly exhibits 
the topological charge 4. This scenario can be explained 
by a change in chirality for a triplet of Dirac points in 
vicinity of the super-merging line (1 — 3-1-3 — 3) — > 
(1 — 3-1-3-1-3). However the behavior of solutions on the 
super-merging line in N8 graphene is beyond the scope 
of the present paper. 

The Nn graphene model was shown to exhibit more 
than one touching point in each valley. Now it remains to 
answer the question whether large Chern number phases 



graphene 


super-merging 
+^ i^, +^ 

^1 C3 ^7 


/(K + q) 


charge 


Nl 


10 


TT* 


1 


N3 


11/2 


TT^ 


-2 


N4 


1 2/5 1/5 


*2 
TT TT 


1 


N7 


1 7/12 1/4 1/12 


7r*7r^ 


-2 



TABLE II. Super-merging characteristics at K. The function 
/ from the effective low-energy Hamiltonian, H^a = |cr+/ -f 
H.c, is written as a function of small momenta tt = + iqy 
and up to a multiplicative constant which is neglected. 

become possible when gapping them with a Haldane 
mass. As long that the position of the band touchings 
and their respective chirality is known, determining the 
topological phase diagram is within analytical grasp. 

IV. CHERN NUMBER PHASE DIAGRAM FOR 
THE LONG DISTANCE HOPPING HALDANE 
MODEL 

The Haldane model is built on the hexagonal lattice 
for Nl graphene by adding N2 (intra-sublattice) hop- 
ping t2, such that when hopping is performed clockwise 
in the unit cell an electron gains a phase cj). However 
there is no net magnetic flux in the unit cell. The N2 
hopping term leads to two contributions of respective 
form ho{k)aQ with ho^k) — /io(— fc) and hj,[k)a^ with 
/i3(fc) = —h^{~k). These two contributions break chiral 
symmetry, but do not break inversion symmetry. The 
first contribution breaks particle-hole symmetry, while 
the second breaks time-reversal symmetry. As far as 
Chern number calculation Eq. ([4| is concerned, the first 
contribution can be discarded and only the second con- 
tribution is relevant, because it produces the necessary 
mass term Eq. ([5| that gaps the Dirac points. As we 
shall see later, when added to Nn graphene model, the 
latter contribution allows for Chern phase diagram with 
only odd (even) Chern number phases. In order to have 
a Chern phase diagram allowing for transition between 
even and odd Chern number phases, it is necessary to 
add a mass term that breaks inversion symmetry. The 
simplest such term is of the form Mas and corresponds 
to a different on-site potential energy on each sublattice. 

The mass term h^as in the original Haldane model 
breaks time-reversal and inversion symmetry. It reads 

/i3 = Af— 2t2 sin(/)[sin(k-a2)— sin(k-ai)+sin(k-(ai— a9))]. 

When intra-sublattice hopping between distant sites is 
allowed, the generalized mass term reads 

/i3 = Af - ^ 2t("^ sin(n0)[sin(7ik • £12) - sin(nk • ai) 

n 

+ sin(nk- (ai -az))], (29) 

where n is an integer that indicates that hopping takes 
place between AA or BB sites situated at a distance of 
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hopping physical distance chemical distance 



t2 

ts 



3 
2^3 



TABLE III. The first hopping integrals tn contributing to the 
Haldane mass. The hopping distances are expressed in units 
of lattice constant. 



n\/3a. Here will be considered only the first two terms 
in this expression, corresponding to a hopping across two 
unit cells (see Fig. [T] and Tab. III). The term containing 
the hopping integral just multiplies the identity Pauli 
matrix and is neglected. Interesting for the topology of 
the problems are hoppings along the links where the elec- 
trons gain the phase (j). Here only the first two terms in 
the mass term are considered, t2 and ig- 

The goal of this part is to illustrate how gapping the 
graphene system with 2n Dirac points can yield Z topo- 
logical phases characterized by a large Chern number (up 
to C = ±n). The following subsections investigate cases 
where different mass term gaps the previously obtained 
Nrt graphene. The strategy will be to illustrate the possi- 
bility of large Chern phases by considering first the action 
of t2 Haldane mass on different models of Nn in Sec. |IV A] 
In Sec. |IVB] it is shown that the addition of terms al- 
lows to further increase the absolute value of the Chern 
number. 



A. t2 Haldane model 

Nl graphene with a hopping t2 constitutes the origi- 
nal Haldane model. The phase diagram is obtained by 
observing that changes sign between the Dirac points 
(±K) of graphene. Therefore the Hamiltonian exhibits 
three topological phases; a trivial insulating phase and 
two C = ±1 QAH phases. Eq. ^ yields in this case 



C = -{sgnM+~sgiiM-). 



(30) 



where M± = M T 3^3*2 sin(/) is the mass term at =fK. 
The phase diagram is represented in Fig. [9j The lines 
M± = represent topological transition lines where the 
bulk gap closes at least at one of the ±K points. 

Larger Chern phases become possible when the under- 
lying model is N3 graphene. Now the mass term takes 
different values between a regular Dirac point and its 
satellites. Therefore the topological charges can add up 
to yield Chern |C| = 2 phases. 

Momentum ztk^ locates any satellite point of ±K and, 
manifestly, the expression for x(k,j) holds in the range of 
existence of separate satellite points. 

Let us define the mass at the regular Dirac points 
M± = M{t (^^'O))- The mass at the satellite Dirac 
points is denoted by m± , if it is associated to the regular 
Dirac points =fK. Then from Eq. ^ it follows that the 




FIG. 9. (Color online). Chern number phase diagram for the 
Haldane Hamiltonian as a function of the on-site energy M 
divided by the hopping integral t2 as a function of the flux 
(f}. The topologically nontrivial insulating phases are color 
identified and have the topological index denoted inside the 
respective regions. The topologically insulating regions, C = 
0, are white. 




FIG. 10. (Color online). Chern number phase diagram for the 
t2 Haldane model on N3 graphene. The hopping parameters 
are t2 — 1/3 and ts = 0.35 in units of ti. 



Chern yields 
1 



C 



(sgnA^_ — sgnA^+) — 3(sgnm__ — sgnm^ 



(31) 



where the mass of the Dirac points read 
M±^MT 373^2 sin 



m± =Mt 2*^(1 +t3)Wl- (^^^ ) sin0. (32) 



Eq. (31 ) yields the phase diagram for the system when 



all 8 Dirac points are present. When there are no satellite 
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FIG. 11. (a) Energy dispersion at the topological transition between C = —2 and C = phases at the merging point between 
the regular — K and its three satellites — ki in N3 graphene. The energy dispersion in N3 Haldane shows a quadratic band 
touching at — K. The parameters are chosen (j> = 7r/2, AI = \/3, ^2 = 1/3 and = 1/2. (b) Energy dispersion for the N3 
Haldane model at the transition between C = 1 and C = —2 phases. The Dirac cones form at the satellite points of — K for 
t2 — 1/3, ts = 0.35 in units of ti. The change in Chern number by three units is reflected in the presence of three Dirac cones 
at the topological transition. 



Dirac points (^3 G (—1, 1/3)), the topology of the system 
is in fact identical to the original system ^3 = and there- 
fore it has the phase diagram in Fig.|9] When ^3 is varied 
to go outside the region (—1, 1/3), two phases of higher 
Chern number develop around the M = line. For exam- 
ple, from Eqs. (32 ), we see that at M = a regular Dirac 



point and its satellites will have the same mass. There- 
fore the Chern number reduces to C = sgnA^+ — sgnA^_ . 
This yields topological phases indexed by ±2. By in- 
creasing |M|, one crosses a transition line where where 
the Haldane mass of all satellite points in the system be- 
comes identical, while it remains different for the regular 
Dirac points. This transition is given by 



m± — 0. 



(33) 



This region extends up to the the last topological transi- 
tion line given by M = ±373^2 sin (/). In this region the 
Chern number reduces again to the original case {tz = 0) 
with C = l/2(sgnA^_ — sgn7W+). When M is increased 
even further, all Dirac points are gapped identically and 
therefore this is the topologically trivial region. In Fig.[TO| 
is represented a typical phase diagram for the case where 
satellite Dirac points are present. 

Note that at the merging point = 1/2 the C = ±1 
phases completely vanish, and the phase C = ±2 would 
have maximal area dehmited by M = ±3%/3t2 sm m. 
Then at the topological transition from |C| = 2 phase 
to the trivial insulator, there is a quadratic band touch- 
ing that is represented in Fig. 11 'a). 



The phase diagram in N3 Haldane model (Fig. 10 1 has 



the nice feature that it accommodates lines of transition 
where Chern number changes by 3 units. This is realized 
by the formation of three Dirac points at the topological 
transition. These band touchings come from the vanish- 
ing of the Haldane mass at the three satellite Dirac points 
previously found in N3 graphene. For example, let us 



take parameters ti = 1, <2 = 1/3 and = 0.35 from the 
phase diagram in Fig. 10 Then fixing (f> — 7r/2, there are 
two transition points between C = —2 and C ~ 1 phases 
near ±K. In particular, near — K, the Dirac points form 
at the satellites where 771+ = 0. The energy dispersion at 
the topological transition is illustrated in Fig. [Tl|b). 

Similarly one can take as underlying model the N4 
graphene model which contains the ti hopping. This was 
shown to produce 7 Dirac points per valley. Hence one 
can expect the presence of larger Chern phases. This is 
exemplified in Fig. |12[ where a choice of particular pa- 
rameters yields |C| = 4 QAH phases. Note the presence 
of multiple Dirac points is reflected in the phase diagram 
as a multiplication of transition lines in M direction for 
fixed magnetic flux cf) (7^ 0, tt). 



B. te Haldane model 

The existence of 2n Dirac points for a submodel con- 
taining only two sigma matrices allows in principle to 
build topological insulators with Chern phases C = n. 
For the N3 graphene model with 8 Dirac points, one 
can have a large Chern number C — ±4. To actualize 
all possible topological phase it is sufficient to add a te 
mass term. It has the effect to produce oscillations in 
the phase dependent Haldane mass, such that the term 
changes sign between a regular graphene Dirac point and 
its satellites in N3 graphene. As expected, all phases are 
attainable under this modification of the Hamiltonian. 

The mass term becomes 

= M — 2<2 sin (/)[sin(k • 3.2) — sin(k • ai) 

-|-sin(k • (ai - a2))] - 2te sin(20)[sin(2k • a2) 
— sin(2k • ai) + sin(2k • (ai — a2))] 

(34) 
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FIG. 13. All QAH phases possible for N3 graphene with 
te Haldane mass; here the phase diagram for the parame- 
ter choice ti — 1, t2 = 1/3, is ~ 0.35, and te — 0.26 illustrates 
this point. For A/ — 0, the possible Chern phases have only 
even Chern number. 



The new phase diagram is computed by considering the 
mass term (34 1 at the 8 N3 graphene Dirac points. Then 



the topological transition hnes are given by the zeroes of 
the new mass terms, A^j_, m^, expressed as a function of 



the previous mass terms from Eq. ( 32 ) 



M'±=M±± 373^6 sin 20 

to'_j_ = m± ^ 2ig sin 20(2 sin 2k — sin 4k) , (35) 

where k — arccos[(t3 — l)/(2t3)] in the domain of exis- 
tence of satelhte Dirac points in N3 graphene. 

The dependence of the mass term on sin 20 makes pos- 
sible large Chern number phases \C\ = ±4 by having 



FIG. 14. (Color online). A Dirac point that is represented by 
• (o) has chirality + (— ). The colored lines represent lines of 
zeros for hi (green), /i2 (red), and the mass term hg (blue). 
The regular Dirac points placed at ( ± , O) are gapped by 
a Haldane mass that has opposite sign. Also the mass term 
changes sign between the regular Dirac point and its satellites. 
For parameters ti = 1, t2 = 1/3, tg = 0.35,i6 = 0.26, M = 
and (j) = n/8 the phase is C = —4. 



the mass term changing sign between the regular Dirac 
cones and its time-reversed one and its own satellites (see 
Fig. 14). When system parameters are varied, N6 Hal- 
dane model can present all Chern phases between —4 and 
4. A phase diagram that illustrates this point is repre- 
sented in Fig. 14 1. The phase diagram was also sampled 



by numerical integration over the BZ in Eq. Q and the 
results were in agreement. 

Let us consider shortly the case of N4 and N7 graphene 
by adding t4, and respectively, hopping terms. With 
hopping integral ti fixed as before, there are two free 
parameters and 14. The parameter space becomes too 
large to describe analytically the dynamics of the Dirac 
points and to track at the same time the sign of the 
mass at the Dirac points. The general thesis remains 
however correct. Larger and larger QAH phases become 
possible. In the case of N4 graphene there is a maximum 
of 6 Dirac points near a K point; for N7 graphene there 
are 9 possible Dirac points per valley. That indicates 
that with a proper mass term one can have the largest 
Chern phases |C| = 7 (in N4 graphene) or |C| = 10 (in N7 
graphene). In Fig. 15 is represented a Haldane mass 
on a N4 graphene with QAH phases |C| < 5. It appears 
that one needs even longer hopping terms in the Haldane 
mass to realize the largest |C| = 7 phase. 

Note that in all the cases the presence of distant- 
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FIG. 15. Haldane model from N4 graphene with a mass 
term. Hopping integrals ti = 1, t2 = 1/3, is = 0.43, t4 — 0.3, 
= 0.35. For M = 0, the possible Chern phases have only 
odd Chern number. 



neighbor hoppings in Haldane mass potentially leads to 
more bulk gap closings at a given on-site energy AI for 
(j) varying from — tt to tt. This is reflected in the struc- 
ture of the phase diagrams, which present oscillations of 
the topological phase boundaries in (j) direction. This 
accounts for the oscillatory nature of the Haldane mass, 
which can pass more times through zero (as a function of 
the flux), when it contains strong distant neighbor hop- 
ping terms. 



V. CONCLUSION 

We have shown that in a graphene-like system adding 
distant-neighbor hopping integrals leads to the appari- 
tion of satellite Dirac points in the spectrum near the 
regular ±K points of graphene. The number of addi- 
tional Dirac points grows as more distant-neighbors are 
considered. Here, Dirac points up to N7 (next x 4 nearest 
neighbor) graphene model are determined. Each new dis- 
tant hopping integral between AB sites potentially yields 
a triplet of Dirac points near K (and because TRI, a 
triplet at — K). For N7 graphene there is a maximum of 
3 triplets of satellites created. 

The position of the nodes in the dispersion requires 
solving a polynomial whose degree grows as more distant- 
neighbors are considered. Analytically, one can hope 
to determine their position only for a limited num- 
ber of added neighbors (here N4 graphene). Already, 
for N4 graphene, the investigation revealed a rich phe- 



nomenology for band touchings in the system. Beside 
Dirac point band touchings, there are semi-Dirac points 
(band touchings with a linear dispersion in one direc- 
tion and quadratic in the other), or higher-energy dis- 
persion points. Among the latter, we show that there is 
a unique super-merging band touching at ±K that can 
be understood from a collision scenario of all possible 
Dirac points under a variation of the hopping integrals. 
Their uniqueness in hopping integral parameter space in- 
dicates that they are extremely unstable. Moreover, the 
peculiarity of this point resides in the fact that is charac- 
terized by a high-energy dispersion, but a low topological 
charge. This is due to the fact that the super-merging 
points result from a union of Dirac points organized in 
triplets with alternating chirality. Numerical and ana- 
lytical investigations revealed also a new phenomena in 
N4 graphene: the formation of Fermi lines for a partic- 
ular choice of parameters. The particular constraints to 
obtain them indicates again that they are unstable band 
touchings. 

The creation of multiple Dirac points is a precondition 
to achieve phases with a large Chern number. This is 
put to test by implementing the Haldane model in the 
distant-neighbor hopping graphene. The Haldane mass 
term gaps the Dirac points such that new QAH phases 
appear. We have presented various Chern number phase 
diagrams to illustrate the role of distant-hoppings in the 
Haldane mass term: the flux dependence allows to re- 
solve neighbor Dirac points with the different chirality 
by gapping them with a opposite mass. Said differently, 
the mass term now changes sign not only between K and 
— K, but also between the satellite created near the reg- 
ular Dirac points. In principle, for 2n Dirac points in the 
modifled graphene, phases with Chern number \C\ = n 
can be created. 

As a flnal remarks concerning these Nn Graphene- 
Haldane models, we stress that we do not claim that such 
long range hopping is relevant to Graphene physics. We 
believe however that the phenomenology of complex band 
touchings and large Chern number phases that appears in 
this two band long range hopping model is rather univer- 
sal and might appear as the effective low energy physics 
of a more realistic nearest-neigbor model with N-orbitals 
or N-atoms per unit cell. In support of this view there 
is a recent worlfSI that establishes a mapping of the low 
energy physics of the bilayer graphene (4-atoms per unit 
cell) with that of N3 graphene near supermerging. 
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